Antiferromagnetic Spin Fluctuations in the Metallic Phase of Quasi-Two-Dimensional 

Organic Superconductors 
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We give a quantitative analysis of the previously published nuclear magnetic resonance (NMR) 
experiments in the k-(ET)2X family of organic charge transfer salts. The temperature dependence 
of the nuclear spin relaxation rate l/Ti, the Knight shift Ks, and the Korringa ratio /C is compared 
to the predictions of the phenomenological spin fluctuation model of Moriya, and Millis, Monien 
and Pines (M-MMP), that has been used extensively to quantify antiferromagnetic spin fluctuations 
in the cuprates. For temperatures above Tnmr — 50 K, the model gives a good quantitative 
description of the data in the metallic phases of several k-(ET)2X materials. These materials 
display antiferromagnetic correlation lengths which increase with decreasing temperature and grow 
to several lattice constants by Tnmr- It is shown that the fact that the dimensionless Korringa 
ratio is much larger than unity is inconsistent with a broad class of theoretical models (such as 
dynamical mean-field theory) which neglects spatial correlations and/or vertex corrections. For 
materials close to the Mott insulating phase the nuclear spin relaxation rate, the Knight shift and 
the Korringa ratio all decrease significantly with decreasing temperature below Tnmr. This cannot 
be described by the M-MMP model and the most natural explanation is that a pseudogap, similar to 
that observed in the underdoped cuprate superconductors, opens up in the density of states below 
Tnmr. Such a pseudogap has recently been predicted to occur in the dimerised organic charge 
transfer salts materials by the resonating valence bond (RVB) theory. We propose specific new 
experiments on organic superconductors to elucidate these issues. For example, measurements to 
see if high magnetic fields or high pressures can be used to close the pseudogap would be extremely 
valuable. 



I. INTRODUCTION 

In the past twenty years a diverse range of new strongly 
correlated electron materials with exotic electronic and 
magnetic properties have been synthesized. Examples in- 
clude high-temperature cuprate superconductors man- 
ganites with colossal magnetoresistance)^ cerium ox- 
ide catalystsr' sodium cobaltates;^ ruthenates^'i^ heavy 
fermion materials^ and superconducting organic charge 
transfer salts. Many of these materials exhibit a sub- 
tle competition between diverse phases: paramagnetic, 
superconducting, insulating, and the different types of 
order associated with charge, spin, orbital, and lattice 
degrees of freedom. These different phases can be ex- 
plored by varying experimental control parameters such 
as temperature, pressure, magnetic field, and chemical 
composition. Although chemically and structurally di- 
verse the properties of these materials are determined 
by some common features; such as, strong interactions 
between the electrons, reduced dimensionality associated 
with a layered crystal structure, large quantum fluctu- 
ations, and competing interactions. Many of these ma- 
terials are characterized by large antiferromagnetic spin 
fluctuations. Nuclear magnetic resonance spectroscopy 
has proven to be a powerful probe of local spin dynamics 
in many strongly correlated electron materialsi^ d°'^^d^ 
The focus of this paper is on understanding what in- 
formation about spin fluctuations can be extracted from 
NMR experiments on the organic charge transfer salts. 

The systems which are the subject of the current 



study are the organic charge transfer salts based on 
electron donor molecules BEDT-TTF (ET), in partic- 
ular the family k-(ET)2X (where k indicates a partic- 
ular polymorph!^). Similar physics occurs in the other 
dimerised polymorphs, such as the /?,/?', and A phases 
These materials display a wide variety of unconventional 
behaviours^ including: antiferromagnetic and spin liquid 
insulating states, unconventional superconductivity, and 
the metallic phase which we focus on in this paper. They 
also share highly anisotropic crystal and band structures. 
However, for various sociological and historical reasons, 
the K salts have been far more extensively studied, and 
because we intend, in this paper, to make detailed com- 
parisons with experimental data, we limit our study to 
K phase salts. This begs the question: do similar phe- 
nomena to those described below occur in the /?, /?', or A 
salts? We would suggest that the answer is probably yes 
but this remains an inviting experimental question. 

The metallic phase of k-{ET)2X is very different from 
a conventional metallic phase. Many features of the 
metallic phase agree well with the predictions of dy- 
namical mean field theory (DMFT)i'* which describes 
the crossover from a 'bad metal' at high temperatures 
to a Fermi liquid as the temperature is loweredj ^^'^^d^ 
This crossover from incoherent to coherent intralayer— 
transport has been observed in a number of experiments 
such as resistivity,^^ thermopower j^^d^ and ultrasonic 
attenuation. The existence of coherent quasiparticles 
is also apparent from the observed magnetic quantum os- 
cillations at low temperatures in k-(ET)2X ^^^i^^^ How- 
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ever, nuclear magnetic resonance experiments (see Figs. 
[T]and[21) on the metallic phase on k-(ET)2X are not con- 
sistent with a Fermi liquid description. The nuclear spin 
relaxation rate per unit temperature, l/TiT, is larger 
than the Korringa form predicted from Fermi liquid the- 
ory. As the temperature is lowered 1 /TiT reaches a maxi- 
mum; we label this temperature Tnmr (the exact value of 
Tnmr varies with the anion X, but typically, Tnmr ~ 50 
K, see Fig. [Ij. l/TiT decreases rapidly as the tempera- 
ture is lowered below Tnmr [see Fig[T] J^i^i^ The Knight 
shift also drops rapidly around rNMR-^ This is clearly in 
contrast to the Korringa-like behavior one would expect 
for a Fermi liquid in which l/TiT and Kg are constant 
for T <C Tp, the Fermi temperature. A similar non- Fermi 
liquid temperature dependence of l/TiT and Ks is ob- 
served in the cuprateSf^MI It has been argued that the 
large enhancement of the measured l/TiT in cuprates 
is associated with the growth of antiferromagnetic spin 
fluctuation within the Cu02 planes as the temperature is 
loweredi ^^'^^ The large decrease observed in l/TiT and 
Kg measurements for underdoped cuprates^ at temper- 
atures well above Tc is suggestive of a depletion of the 
density of states (DOS) at the Fermi level which might 
be expected if a pseudogap opens at Tnmr- 

A quantitative description of spin fluctuations in the 
metallic phase of k-(ET)2X has not been given pre- 
viously. However, the importance of spin fluctuations 
for the superconducting k-(ET)2X has been pointed out 
by several groups . ^'^"I'^^'^^i'^^'^'^i'^^ Since superconductiv- 
ity arises from an instability of the metallic phase, it is 
important to understand the strength of the spin fluctu- 
ations in the metallic phase. 

We use the phenomenological antiferromagnetic spin 
fluctuation model which was first introduced by Moriya 
in his self consistent renormalization (SCR) theory^^ and 
then appHed by Millis, Monien and Pines (MMP)^^ to 
cuprates, to examine the role of spin fluctuations in the 
metallic phase of k-{KT)2X . We fit the spin fluctuation 
model to the nuclear spin relaxation rate per unit temper- 
ature 1/TiT, Knight shift Kg , and Korringa ratio JC data. 
We find that the large enhancements measured in 1/TiT 
and IC above Tnmr are the result of large antiferromag- 
netic spin fluctuations [see Figs. [T]and[2]. The antifer- 
romagnetic correlation length increases as temperature 
decreases and the relevant correlation length is found to 
be 2.8 ± 1.8 lattice spacings at T = Tnmr = 50 K. The 
model produces a reasonable agreement with experimen- 
tal data down to T '--^ 50 K. The spin fluctuation model 
predicts a monotonically increasing 1/TiT with decreas- 
ing temperature while the measured 1/TiT below 50 K 
is suppressed but never saturates to a constant value. 
This is contrary to what is expected for a Fermi liquid 
where 1/TiT is constant. This indicates that the metallic 
phase of K-(ET)2Ar is richer than a renormalized Fermi 
liquid as has been previously thought to describe the low 
temperature metallic state. 

The structure of the paper is as follows. In Section II 
we introduce the temperature dependence of the nuclear 



spin relaxation rate. Knight shift, and Korringa ratio and 
describe how they probe the dynamic susceptibility. We 
calculate these properties in a number of approximations 
and contrast the results. In Section III we demonstrate 
that the spin fluctuation model provides reasonable fits 
to the existing experimental results for k-{FiT)2X above 
Tnmr and discuss its limitations when applied to those 
materials. In Section IV we discuss the unresolved issues 
and suggest new experiments to understand those issues. 
Finally, we give our conclusions in Section V. 



II. THE SPIN LATTICE RELAXATION RATE, 
KNIGHT SHIFT, AND KORRINGA RATIO 

In this section we discuss the temperature dependence 
of the nuclear spin lattice relaxation rate 1 /Ti , Knight 
shift Kg, Korringa ratio /C, and their dependence on the 
dynamic susceptibility x(q, w) = x'(q, w) -I- ix"(q, w). 
The general expressions for 1/Ti, Kg and JC are given 
by^ 



and 
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IC = 



iim^5:iA(q)i^^::iHii^, (la) 



.0 72^,4 

l^(o)lx'(o,o) ^ 



(lb) 



(Ic) 



where A{q) is the hyperfine coupling between the nuclear 
and electron spins, and (7e) is the nuclear (electronic) 
gyromagnetic ratio. For simplicity we will consider a mo- 
mentum independent hyperfine coupling |A| in what fol- 
lows. Note that Eqs. ^ show that this is an approxima- 
tion for Ti but that it is not an approximation at all for 
Kg ■ This is because Kg only probes the long wavelength 
physics and hence only depends on A{0), the hyperfine 
coupfing at q = 0. 

The calculation of the quantities in Eqs. ([T]) boils down 
to determining the appropriate form of the dynamic sus- 
ceptibility. Below we discuss, in some detail, the dynamic 
susceptibility within the spin fiuctuation model and cal- 
culate 1/TiT, Kg, and IC. The results from dynamical 
mean field theory (DMFT) will also be discussed for com- 
parison. 



A. The Spin Fluctuation Model 

The dynamic susceptibility in this model is given 

by28i29 



X(q, ^) = XLw(w) + XAF(q, w), 



(2) 



where Xhwi^) is the dynamic susceptibility in the long 
wavelength regime and XAF(q, i-^) is a contribution to the 
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dynamic susceptibility which is peaked at some wave vec- 
tor Q. These susceptibihties take the form 



XAF(q,t^) = 



xo(r) 



1 -iuj/T{T) 



xq{t) 



l+^(r)2|q-Q|2-zc./c.sF(r)^ 



(3) 



where Xo{T) [xq{T)\ is the static spin susceptibihty at 
q = [Q], r(r) [a;sF(r)] is the characteristic spin fluc- 
tuation energy which represents damping in the system 
near q = [Q], and ^(T) is the temperature dependent 
correlation length. Hence, the real and imaginary parts 
of the dynamic susceptibility can then be written as 



X'(q,0) 

x"(q,^) 



Xo(T) 

^xo(r) 



^ I xq{t) 



1 



xq{t)t 



1 



xo(r)c.sF(r) (i-f5(r)2|q-Q|2)2 



(4) 

Note that the above form of xlw{'^) is the appropri- 
ate form for a Fermi liquid. Therefore, if the system 
under discussion is not a Fermi liquid then the valid- 
ity of this expression for xlw{^) cannot be guaran- 
teed. For example, the marginal Fermi liquid theory 
predicts a different frequency dependence. If the dy- 
namic susceptibility has a large peak at Q 7^ then 
1/Ti will not be strongly dependent on the long wave- 
length physics [because 1/Ti measures the susceptibility 
over the entire Brillouin zone, c.f., Eq. (fTa|) . and there- 
fore will be dominated by the physics at q = Q]. On 
the other hand, the Knight shift is a measure of the long 
wavelength properties [c.f., Eq. (jlbp ] and therefore may 
be sensitive to the details of xim{'^)- Below we follow 
MMB^ and explicitly assume that the uniform suscep- 
tibility (xo) a-iid the spin fluctuation energy near q = 
(r) are temperature independent. One justification for 
this approximation in organics is that the Knight shift is 
not strongly temperature dependent;^ However, this ap- 
proximation breaks down in systems where the uniform 
susceptibility is strongly temperature dependent such as 
YBaaCusOe.ea^ and Lai.8Sro.i5Cu04^ 

In the critical region ^(T) ^ a, where a is the lattice 
constant, one has^^ 



Xq{T) = 

WSF(r) = 



Xo 



(5) 



where 77 is the critical exponent which governs the power- 
law decay of the spin correlation function at the critical 
point, z is the dynamical critical exponent, and is 
a temperature independent length scale. The simplest 



assumptions are relaxation dynamics for the spin fluctu- 
ations (characterized by z = 2) and mean field scaling 
of the spin correlations [rj — 0). Within these approx- 
imations the real and imaginary parts of the dynamic 
susceptibility are given by 



x'(q,o) 

X"(q,^) 



Xo 
wxo 

r 



1 



[l + e(T)2|q-Q|^ 



l-f /5 



[l + e(r)2|q-Q 



212 



(6) 



where (3 = (ct/^o)^- The temperature independent, di- 
mensionless parameter (3 can also be expressed in terms 
of the original variables appearing in the dynamic sus- 
ceptibility in Eq. ([3]) as 



Xq{T)V ( a 



Xot^SF(T) \£,{T) 



(7) 



Written in this form, [3 has a clear interpretation: it rep- 
resents the strength of the spin fluctuations at the wave 
vector Q relative to those at q 0. We will now consider 
two cases: antiferromagnetic and ferromagnetic spin fluc- 
tuations. 



1. Antiferromagnetic Spin Fluctuations 

If we have antiferromagnetic spin fluctuations then the 
dynamic susceptibility x(q, w) is peaked at a finite wave 
vector q = Q; for example, on a square lattice with near- 
est neighbor exchange only, Q — (7r,7r). The NMR re- 
laxation rate. Knight shift, and Korringa ratio can be 
calculated straightforwardly from the real and imaginary 
parts of the dynamic susceptibility given in Eq. ^ . The 
results are 



1 



l^lxo 



/C = 



1 + mm? 



f3 



1 + [Qm? 



2rxo 



(8a) 
(8b) 

(8c) 



where Q is a cutoff from the momentum integration 
[c.f. Eq. (HH)]. For C{T) > a: 1/TiT ~ ^{Tf, 
and Kg ~ constant which leads to the Korringa ra- 
tio /C ~ (?i72/2rxo)[Q^(T)]2. In this model the Kor- 
ringa ratio can only be equal to unity if the spin fluc- 
tuations are completely suppressed (/3 = 0). Hence, 
one expects /C > 1 if antiferromagnetic fluctuations are 
dominant It has been shown^^, that the Korringa 
ratio is unity when the hyperfine coupling A{q) is mo- 
mentum independent and the vertex corrections are neg- 
ligible. The fact that the Korringa ratio is larger than 
one indicates that there are significant vertex corrections 
when there are large antiferromagnetic fiuctuations. 
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2. Ferromagnetic Spin Fluctuations 

For ferromagnetic spin fluctuations, x(q, cj) is peaked 
at q = 0. The NMR relaxation rate is exactly the same 
as that given in Eq. l|5a|) because 1/TiT comes from 
summing the contributions form all wave vectors in the 
first Brillouin zone, which makes the location of the peak 
in x(qj in the momentum space irrelevant. In contrast, 
the Knight shift will be different in the ferromagnetic and 
antiferromagnetic cases because Ks only measures the 
q = part of the dynamic susceptibility; hence Ks will 
be enhanced by the ferromagnetic fluctuations. Thus, 
for ferromagnetic spin fluctuation description the Knight 
shift Ks is given by 



Ks 



\A\xo 



1 + 



and the corresponding Korringa ratio by 



nil 



1 + /3 



K(r)/a]- 



i+[QC(T)]^ 



2rxo [i + V?(^/a)2 



(9) 



(10) 



For ^(T) > a: l/T{r - £,{Tf, and Ks - ^(r)^ which 
leads to /C ~ {h-il/2Txo)[<{T)/a]-'^. Thus we see that 
/C < 1 in the presence of ferromagnetic fluctuations42ii^ 
So again vertex corrections are important if the system 
has strong ferromagnetic fluctuations. Recall that, in 
contrast, for antiferromagnetic fluctuations the Korringa 
ratio is larger than one. Thus evaluating the Korringa 
ratio allows one to determine whether antiferromagnetic 
or ferromagnetic spin fluctuations are dominant. 



B. Dynamical Mean Field Theory 

DMFT is an approach based on a mapping of the 
Hubbard model onto a self-consistently embedded An- 
derson impurity model»iii^i^ DMFT predicts that the 
metallic phase of the Hubbard model has two regimes 
with a crossover from one to the other at a tempera- 
ture Tq. For T<To the system is a renormalized Fermi 
liquid characterized by Korringa-like temperature de- 
pendence of 1/TiT and coherent intralayer transport. 
Above To, the system exhibits anomalous properties with 
l/TiT ~ a+b{To/T) (c.f., Ref. H) and incoherent charge 
transport. This regime is often refereed to as the 'bad 
metal' Microscopically the bad metal is characterized 
by quasi-localized electrons and the absence of quasipar- 
ticles. This temperature dependence is similar to that 
for the single impurity Anderson niodel4^ Note that this 
temperature dependence is similar to that found for spin 
fluctuations [c.f., Eq. p^ ]. 

The predictions of DMFT correctly describe the prop- 
erties of a range of transport and thermodynamic ex- 
periments on the organic charge transfer salts i^i^^ii^iiii^ 
This suggests that these systems undergo a crossover 
from a bad metal regime for T>Tq to a renormalized 



Fermi liquid below Tq. However, we will show below (also 
see Fig. [1]) that the nuclear spin relaxation rate is sup- 
pressed but never saturates below Tnmr; this is not cap- 
tured by DMFT. This suggests that the low-temperature 
regime of K-(ET)2Ar is more complicated than the renor- 
malized Fermi liquid predicted by DMFT which, until 
now, has been widely believed to be the correct descrip- 
tion of the low temperature metallic state in the organic 
charge transfer salts. 



III. SPIN FLUCTUATIONS IN «:-(ET)2X 

The NMR relaxation rate. Knight shift, and Korringa 
ratio in the antiferromagnetic spin fluctuations model 
are given by Eqs. (|5a)l. (|8b)). and ((5c)l . Their temper- 
ature dependence comes through the antiferromagnetic 
correlation length. We adopt the form of ^(T) from M- 
MMB28i^: C(r)/C(T^) = y/2T^/{T + T^). For this form 
of the correlation length, represents a characteristic 
temperature scale of the spin fluctuations and ^(T) is 
only weakly temperature dependent for T <C TJc. For 
this choice of ^(T) we have 



(TiT)o 
TiT 

Ks 



1 + 

{Ks)o 
1 

/Co- 



(T/T, + 1)2 + 2Tr^C{T/n + 1) 



1 + 27r2C + T/T^ 



(T/T^ + l)^+2Tr^C(T/T^ + l) 



1+2tt^C+T/T^ 



(11a) 
(lib) 

(He) 



where we have defined 



C 

(i/riT)o 

(^.)o 
and /Co 

to simplify the notation. 



2TrkB\A\^Xo 



2rxo' 



(12) 



A. The Nuclear Spin Relaxation Rate 

We now analyze the temperature dependence of 1/Ti. 
In the discussion to follow, we will assume that the corre- 
lation length is sufficiently large compared to the lattice 
spacing and that the quantity 27r2C = 47r2[^(Ta;)/a]^ is 
much larger than T/T^- These two assumptions imply 
that 27r2C(r/Tr + 1) is more dominant than {T/T^ + if 
in the denominator of the second term inside the square 
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bracket of Eq. (jllap . Keeping only the dominant term, 
we arrive at the expression for 1/TiT 



(riT)o 



1 



1 



T/T, + 1 



(13) 



The assumption 2'k'^C{T/T^ + 1) > (T/T^ + if can be 
easily worked out to give a self consistency condition 



> 



T/T, + 1 
47r2 ■ 



(14) 



We will use this relation later in Section IIIB as one of 
the tests for the validity of our approximation. 

The NMR relaxation rate per unit temperature calcu- 
lated from the spin fluctuation model [c.f., Eq. ([8|)] is a 
monotonic decreasing function of temperature. Thus one 
realizes immediately that the data, reproduced in Fig. [l] 
for temperatures below Tnmr is not consistent with the 
predictions of the spin fluctuation theory. We will return 
to discuss this regime latter. We begin by investigating 
the high temperature regime, T > Tnmr- 

We flt the 1/TiT expression, Eq. to the experi- 

mental data of De SotoH^ for K-(ET)2Cu[N(CN)2]Br be- 
tween Tnmr and 300 K with (l/Tir)o, (3[^{T^) / a]^ , and 
Tj. as free parameters. It is not possible to obtain /3 
and ^{Tx)/a independently from fltting to l/TiT data 
because the model depends sensitively only on the prod- 
uct l3[^{Tx)/a]'^ (see Eq. (HH)). The results are plotted 
in Fig. [T] and the parameters from the fits are tabu- 
lated in Table HI We have checked the validity of our 
approximation by plotting \/TiT given by Eq. (jllap for 
K-(ET)2Cu[N(CN)2]Br in Fig. [it, where there is Kor- 
ringa ratio data (see Fig. \^ and thus we can determine 
/3 and ^{Tx)/a individually. It can be seen from Fig. [T]3 
that the disagreement between 1/TiT plotted from Eqs. 
(|llap and (fT3| is smaller than the thickness of the curves. 
Therefore, this approximation is well justified. It will also 
be shown in Section IIIB that the correlation length is in- 
deed rather large and the self consistency condition, Eq. 
p4p . is satisfied, thus providing further justification for 
the use of Eq. (fT^ here. 

The model produces a reasonably good fit to the ex- 
perimental data on K-(ET)2Cu[N(CN)2]Br (Ref. SI) be- 
tween Tnmr, the temperature at which 1/TiT is maxi- 
mum, and room temperature. In the high temperature 
regime (e.g., around room temperature), 1/TiT has a 
very weak temperature dependence, indicating weakly 
correlated spins. The large enhancement of 1/TiT can 
be understood in terms of the growth of the spin fluctu- 
ations: as the system cools down, the spin-spin correla- 
tions grow stronger which allows the nuclear spins to re- 
lax faster by transferring energy to the rest of the spin de- 
grees of freedom via these spin fluctuations. Strong spin 
fluctuations, measured by large values of f3[^{Tx) / a]'^ , 
are not only present in K-(ET)2Cu[N(CN)2]Br but also 
observed in other materials such as fully deuterated k- 
(ET)2Cu[N(CN)2]Br {which wiU be denoted by K(d8)- 
(ET)2Cu[N(CN)2]Br} and k-(ET)2Cu(NCS)2. The re- 



sults of the fits for K(d8)-(ET)2Cu[N(CN)2]Br and k- 
(ET)2Cu(NCS)2 are shown in Fig. [TJ The parameters 
that produce the best fits are also tabulated in Table [J 
In all of the cases studied here, strong spin fluctuations 
are evident from the large value of l3[£,{Tx) / 0]"^ . 

The nature of the spin fluctuations, i.e., whether they 
are antiferromagnetic or ferromagnetic, cannot, even in 
principle, be determined from the analysis on l/TiT. 
Both cases yield the same 1/TiT [see Eq. ((8a)) and Sec 
II. B. 2] because the nuclear spin relaxation rate is ob- 
tained by summing all wave vector contribution in the 
first Brillouin zone. However, in the next section we will 
use the Korringa ratio to show that the spin fluctuations 
are antiferromagnetic. 

Below Tnmr, the calculated 1/TiT continues to rise 
while the experimental data show a decrease in the nu- 
clear spin relaxation rate per unit temperature. However, 
the data do not reach a constant 1/TiT as expected for 
a Fermi liquid. This indicates that the physics below 
Tnmr is dominated by some other mechanism not cap- 
tured by the spin fluctuation theory, Fermi liquid theory, 
or DMFT. 

One might argue that the discrepancy between the the- 
ory and experiments below Tnmr stems from our assump- 
tion of a q- independent hyperfine coupling in the 1/TiT 
expression. However, in section IIII Dl we will show that 
the Knight shift is also inconsistent with the predictions 
of the spin fluctuation model below Tnmr- While in- 
cluding the appropriate q-dependent hyperflne coupling 
might change the temperature dependence of 1/TiT, it 
certainly cannot affect the temperature dependence of 
the Knight shift because Kg only depends on A{0) [as 
can be seen from Eq. (fTb]) ]. 



B. The Korringa Ratio 

In the previous section we compared the predictions of 
the spin fluctuation model for 1/TiT to the experimental 
data and obtained good agreement with the data between 
Tnmr and 300 K. However, we were not able to determine 
(3 and ^{Tx)/a independently because 1/TiT is sensitive 
only to the product P[£,{Tx) / a]'^ . We were also unable 
to determine whether antiferromagnetic or ferromagnetic 
spin fluctuations are dominant. We resolve these ques- 
tions by studying the Korringa ratio JC. It has previously 
been pointed out that antiferromagnetic (ferromagnetic) 
fluctuations produce a Korringa ratio that is larger (less) 
than one.^^J^ We have also seen in Section II that in the 
limit of large correlation lengths, JC ~ (S./a)'^ > 1 for 
antiferromagnetic spin fluctuations and /C ~ (a/C)^ < 1 
for ferromagnetic spin fluctuations. The Korringa ratio 
data for K-(ET)2Cu[N(CN)2]Br (see Fig O arc signifl- 
cantly larger than one at all temperatures which shows 
that antiferromagnetic fluctuations dominate. With this 
in mind, we study the antiferromagnetic spin fluctuation 
model. 

First we note that K^, given by Eq. (jllbp . has a 
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50 100 150 200 250 300 50 100 150 200 250 300 

T(K) T(K) 

FIG. 1: [Color online] Comparison of the measured nuclear spin relaxation rate per unit temperature, 1/T\T, with the predic- 
tions of the spin fluctuation model for various organic charge transfer salts. Panel (a) shows data for K-(ET)2Cu[N(CN)2]Br mea- 
sured by Mayaffre et al.^^. Panel (b) shows data for K-(ET)2Cu[N(CN)2]Br measured by De Soto et al^. Panel (c) shows data 
for fi:(d8)-(ET)2Cu[N(CN)2]Br measured by Miyagawa et al^. Panel (d) shows data for a k-(ET)2Cu(NCS)2 powder sample 
measured by Kawamoto et al}^ The 1/TiT data are weakly temperature dependent at high temperatures, have a maximum 
at Tnmr ~ 50 K, and drop abruptly below Tnmr ~ 50 K, contrary to what one would expect for a Fermi liquid in which 
1/TiT is constant. The remarkable similarities of these data result from the quantitative and qualitative similarity of the 
antiferromagnetic spin fluctuations in the metallic phases of these materials. The parameters that produce the best fits (solid 
lines) to Eq. (|llc|) are tabulated in Table |T] The spin fluctuation model gives a good flt to the experimental data between 
Tnmr ~ 50 and room temperature which suggests strong spin fluctuations in the metallic states of K-(ET)2Cu[N(CN)2]Br, 
K(d8)-(ET)2Cu[N(CN)2]Br, and k-(ET)2Cu(NCS)2. However, below Tnmr the spin fluctuation model does not describe the 
data well, indicating that some other physics dominates over the spin fluctuations physics. In each figure the solid line is 
obtained from the approximate form for 1/TiT given by Eq. (|13p . To check that this approximation is reasonable, we also plot 
1/TiT without any approximation, given by Eq. (|llap . as a dashed line in panel (b). The full and dashed lines cannot be dis- 
tinguished until well below Tnmr and so we concluded that the approximation is excellent in the relevant regime. Note that the 
analysis on 1/TiT cannot differentiate between antiferromagnetic and ferromagnetic spin fluctuations (see section Hi A 2|l . but 
the Korringa ratio strongly differentiates between these two case and indicates that the fluctuations are antiferromagnetic (see 
Fig. [2)1. The nomenclature K-Br, d8-Br, and k-NCS is used as shorthand for K-(ET)2Cu[N(CN)2]Br, K(d8)-(ET)2Cu[N(CN)2]Br, 
and K-(ET)2Cu(NCS)2 respectively in the figure keys. 



weak temperature dependence because of our assump- 
tion that [^{Tx)/a] is generally larger than unity and 
27r^C ^ T/Tx- Thus, the second term inside the 
square bracket in Eq. (jllbp can be approximated by 
v^C/(H-27r2C + T/T^)-i ~ \/3/(27r2) and the Knight 
shift wiU be given by A% ~ {Ks)o[l + v^/(27r2)] which 
is temperature independent. We use this temperature 
independent Knight shift to calculate the Korringa ratio 



ICn 1 



7e y 1 



(15) 



7r2(T/T, + l); Vl + \/?/(27r2) 

where the prefactor /Co is given by Eq. 

We fit Eq. to the experimental data for the Kor- 
ringa for K-(ET)2Cu[N(CN)2]Br.25 The result is plotted 
in Fig. O The Korringa ratio data are well reproduced 
by the antiferromagnetic spin fluctuation model when 
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Material 


Ref. 


(i/rir)o (s-^K-i) 


(K) 


Tnmr (K) 




^(Tnmr)/^ 


K-Bl 


Mayaffre [24] 


0.09 ± 0.01 


7 ±6 


60 


290 ±250 


7 


K-Br 


De Soto [25] 


0.02 ± 0.01 


20 ±10 


50 


680 ±430 


2.8 ± 1.8 


d8-Br 


Miyagawa [471 


0.04 ±0.01 


6 ±4 


40 


85 ±65 


? 


K-NCS 


Kawamoto [48] 


0.06 ±0.01 


11 ±3 


55 


110 ±90 


? 



TABLE I: The parameters obtained from the fits which are used to produce Fig. [T] Evidence for strong spin fluctuations 
come from the large value of P[£,{Tx)/a]'^ which are present for all the materials tabulated above. Tnmr is determined from 
the peak of 1/TiT [see Fig. [T]. The correlation length shown in the last column in the table was obtained by analyzing 
the Korringa ratio data available for K-(ET)2Cu[N(CN)2]Br [see Section IIIB]. The correlation length for Mayaffre K-Br,— 
Miyagawa dS-Br,'*'' and Kawamoto's k-NCS^^ can not be determined from our analysis because there are not sufficient data. 
This is shown with question marks in the correlation length column. In the table K-Br, d8-Br, and k-NCS are used as shorthand 
for K-(ET)2Cu[N(CN)2]Br, K(d8)-(ET)2Cu[N(CN)2]Br, and k-(ET)2Cu(NCS)2 respectively. 




2 I , , , , , ] 

50 100 150 200 250 300 
T(K) 



FIG. 2: [Color online] Comparison of the Korringa ratio 
K. <x l/TiTKl of K-(ET)2Cu[N(CN)2]Br measured by De 
Soto et al?^ with the prediction of the antiferromagnetic spin 
fluctuation model. The best fit to Eq. (|15p is indicated by the 
solid line. The Korringa ratio is larger than 1 which indicates 
that the spin fluctuations are antiferromagnetic (/C < 1 for 
ferromagnetic fluctuations, see Section III A 2\ . The antifer- 
romagnetic correlation length is found to be 2.8 ± 1.8 lattice 
spacings at T = 50 K. Below T = 50 K the Korringa ratio is 
suppressed, and the spin fluctuation model does not explain 
this behavior. This is a clear indication that different physics 
is at play below 50 K. 



T > Tnmr- This is again consistent with our earlier con- 
clusion that the spin fluctuations are antiferromagnetic. 
In this fit we have three free parameters, /3[^(Ta;)/a]^, T^, 
and \/^, two of which, [3[^{Tx) / and T^, have been de- 
termined from fitting 1/TiT. There is only one remaining 
free parameter in the model, which can then be de- 
termined unambiguously from the Korringa fit yielding 
/3 = 60 ± 20. This value of P implies that the antiferro- 
magnetic correlation length £,{T) — 2.8 ± 1.8a (a is the 
unit of one lattice constant) at T = 50 K. This value 
is in the same order of magnitude as the value of the 
correlation length estimated in the cupratesi^ 

We now return to discuss the validity of our approx- 



imation which was stated in the beginning of Section 
III A. The correlation length has been determined to be 
^{T)/a = 3.7 ± 2.4 at T = = 20 K from the fit to the 
Korringa ratio data. This result surely satisfies the re- 
quirement that the correlation length is larger than unit 
lattice spacing. A stronger justification for our approxi- 
mation comes from the self consistency relation Eq. (jl4p . 
With £,{T^)/a = 3.7 and = 20 K, one could easily 
check that Eq. ((T4)) is indeed satisfied in the relevant 
regime, i.e. between Tnmr and room temperature. 

A large Korringa ratio^i^ has previously been ob- 
served in the cuprates indicating similar antiferromag- 
netic fluctuations in these systems. The Korringa ratio 
has also been measured in a number of heavy fermion 
compounds j^^'^^1^^ Similar antiferromagnetic fluctua- 
tions are also present in CeCu2Si2; the Korringa ratio of 
this material has a value of 4.6 at T = 100 mK (Ref. 's?). 
In contrast, YbRh2Si2^ and CeRu2Si2^'^, show strong 
ferromagnetic spin fluctuations as is evident from the 
Korringa ratio less than unity. In Sr2Ru04^ the Ko- 
rringa ratio is approximately 1.5 at T = 1.4 K. Upon 
doping with Ca to form Sr2_a;Caa;2Ru04, the Korringa 
ratio becomes less than one which indicates that there 
is a subtle competition between antiferromagnetic and 
ferromagnetic fluctuations in these ruthenates. 



C. The Antiferromagnetic Correlation Length 

It is important to realize that the spin fluctuation for- 
malism can be used to extract quantitative information 
about the spin correlations from NMR data. From the 
fit for K-(ET)2Cu[N(CN)2]Br (Table|I| we found that the 
antiferromagnetic correlation length ^(T)/a = 2.8 ±1.8 
at T = 50 K. In order to understand the physical signif- 
icance of this value of ^ (T) it is informative to compare 
this value with the correlation length for the square^^ and 
triangular^® lattice antiferromagnetic Heisenberg models 
with nearest neighbor interaction only. 

It has been shown^ that, on the square lattice, the 
antiferromagnetic Heisenberg model with nearest neigh- 
bor interaction only has a correlation length of order 
$(T)/a ~ 1 for T = J and of order i{T)/a - 30 for 
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T = 0.3J. On the other hand for the antiferromag- 
nctic Heisenberg model with nearest neighbor interac- 
tion only on the isotropic triangular lattice, the corre- 
lation length is only of order a lattice constant at T = 
0.3J.5'5 Thus the correlation length, ^iT)/a = 2.8 ± 1.8 
at T — 50 K, obtained from the analysis of the data 
for K-(ET)2Cu[N(CN)2]Br is reasonable and places the 
materials between the square lattice and isotropic tri- 
angular lattice antiferromagnetic Heisenberg model as 
has been argued on the basis of electronic structure 
calculations 

One of the best ways to measure antiferromagnetic 
correlation length is by inelastic neutron scattering ex- 
periments. To perform this experiment, one needs high 
quality single crystals. Unfortunately, it is difhcult to 
grow sufhciently large single crystals for k-(ET)2X; how- 
ever, recently some significant progress has been made in 
this direction.^* Another way to probe the correlation 
length is through the spin echo experiment. The spin 
echo decay rate I/T2 is proportional to the temperature 
dependence correlation length. To the authors' knowl- 
edge there is no spin echo decay rate measurement on 
the metallic phase of the layered organic materials at the 
present time. Thus, it is very desirable to have such ex- 
perimental data to compare with the value of ^{T) we 
have extracted above. 



D. The Knight Shift 

As we pointed out in Section II the Knight shift Kg will 
generally have a weak temperature dependence through- 
out the whole temperature range and so, thus far, we 
have neglected its temperature dependence. However, it 
is apparent from Eq. (jllcp that for any choice of param- 
eter values {P, ^{Tx)/a, and T^}, Kg will always increase 
monotically as the temperature decreases. Therefore the 
temperature dependence of the Knight shift potentially 
provides an important check on the validity of the spin 
fluctuation model. However, in the following discussion 
one should recall the caveats (discussed in section III A[) 
on the validity of the calculation of the Knight shift stem- 
ming from the assumption that the dynamics of the long 
wavelength part of dynamical susceptibility relax in the 
same manner as a Fermi liquid does. 

In contrast to the prediction of the spin fluc- 
tuation model, the experimental data^^ for k- 
(ET)2Cu[N(CN)2]Br show that K,, decreases slowly 
with decreasing temperature which then undergoes a 
large suppression around Tk^ ~ 50 K. It should be 
emphasized here that Tk^ is approximately the same as 
Tnmr, the temperature at which 1/TiT is maximum. 

Since it is not possible to explain any of the NMR data 
below Tnmr in terms of the spin fluctuation model within 
the approximations discussed thus far, we focus on the 
temperature range between 50 K to 300 K just as we 
did for the analysis of l/TiT. Even in this temperature 
range, there is a puzzling discrepancy between theory and 



experiment: the experimental data decrease slowly with 
decreasing temperature while the theoretical calculation 
predicts the opposite. We will argue below that this dis- 
crepancy arises because the data are obtained at con- 
stant pressure while the theoretical prediction assumes 
constant volume. Since the organic charge transfer salts 
are particularly soft, thermal expansion of the unit cell 
may produce a sizeable effect to the Knight shift and may 
not be neglected. In principle, an estimate of the size of 
this effect could be made following Wzietek et al.^ as 

^Ks = Kl-Kl (16) 

where is the (experimentally obtained) isobaric 
Knight shift, is the (calculated) constant volume 
Knight shift, {VdP/dV)T is the isothermal compress- 
ibility, and ^{dV/dT)p is the linear thermal expansion. 
However, it is not possible to obtain an accurate esti- 
mate for Aifs at this time because there are no complete 
data sets for iff, isothermal compressibility, and ther- 
mal expansion as a function of temperature and pressure 
for the K-(ET)2-X family. However, a rough estimate 
for AA's may be made using the available experimental 
data4^ This suggests that the experimental data is con- 
sistent with the spin fluctuation theory. Clearly, further 
experiments are required to test this claim conclusively. 
Therefore we raise this issue predominately to stress the 
importance of systematic measurements of the parame- 
ters in Eq. (fTB]) . 

Given the large uncertainty in AKs we take Kg to 
be constant for temperatures above 50 K in the rest of 
this paper. This is clearly the simplest assumption, it 
is not (yet) contradicted by experimental data, and, per- 
haps most important, any temperature dependence in the 
Knight shift is signiflcantly smaller than the temperature 
dependence of 1/TiT. 

Regardless of the size of AXs , the Knight shift calcu- 
lated from the spin fluctuation model is inconsistent with 
the experimental data below Tr-^ ~ 50 K (see Fig. 4 in 
Ref. Wa ). The calculated Kg shows a weakly increasing 
Ks with decreasing temperature, while the measured 
is heavily suppressed below 50 K. One important point 
to emphasize here is that the temperature dependence of 
Ks will not change even if one uses the fully q-dependent 
v4(q) since Ks only probes the q = component of 
the hyperfine coupHng and susceptibility [see Eq. (fTbl) ]. 
Thus, putting an appropriate q-dependent hyperfine cou- 
pling will not change the result for Kg (although it might 
give a better description for l/TiT). This provides a 
compelling clue that some non-trivial mechanism is re- 
sponsible to the suppression of l/TiT, Ks, and K. below 
50 K. 

We have not addressed how the nuclear spin relaxation 
rate is modified by the thermal expansion of the lattice. 
Since the organic compound is soft, it is interesting to 
ask if there is a sizeable effect to 1/TiT. Wzietek et 
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al^ have performed this analysis on quasi-lD organic 
compounds whose relaxation rate in found to scale like 
Xg. One can straightforwardly derive the effect of vol- 
ume changes from the Hubbard model. If one uses the 
relation 1/TiT ^ Xs ^'^d assumes fixed U and i, then 
l/TiT - 1/^2 wiU follow. However, it is clear fr om the 
phase diagram of the organic charge transfer salts [see 
Ref. 8| that there is a rather large change in U and t 
for even small pressure variations. Therefore, there is 
no obvious relationship between 1/TiT and Xs for the 
quasi-2D organics and it is not clear how the imaginary 
part of the susceptibility x"(q, a;), which enters l/TiT, 
is effected by thermal expansion and lattice isothermal 
compressibility. Again, this stresses the importance of 
the detailed experiments needed to determine the effect 
of thermal expansion of the lattice on the measured re- 
laxation rate. 



IV. OPEN PROBLEMS AND FUTURE 
EXPERIMENTS 

Open problems. The large suppression of l/TiT and 
Ks below Tnmr observed in all the k salts studied here 
cannot be explained by the M-MMP spin fluctuation 
model. One plausible mechanism to account for this fea- 
ture is the appearance of a pseudogap which causes the 
suppression of the density of states at the Fermi energy. 
This is because at low temperature l/TiT and Kg are 
proportional to fP{Ep) and p{Ep), where p{Ep) is the 
full interacting density of states at the Fermi energy42i 
Independent evidence for the suppression of density of 
states at the Fermi level comes from the linear coeffi- 
cient of the specific heat The electronic specific heat 
probes the density of excitations within fc^T of the Fermi 
energy. Any gap will suppress the density of states near 
the Fermi surface which results in the depression of the 
specific heat coefficient 7. Kanoda^ compared 7 for sev- 
eral K-(ET)2A salts and found that in the region close 
to the Mott transition, 7 is indeed reduced. One possi- 
ble interpretation of this behavior is a pseudogap which 
becomes bigger as one approaches the Mott transition. 
However, other interpretations are also possible. In par- 
ticular one needs to take care to account for the possible 
coexistence of metallic and insulating phases; this is ex- 
pected as the Mott transition is first order in the organic 
charge transfer salts. ^^'^^ The existence of a pseudogap 
has also been suggested in A-(BEDT-TSF)2GaCl4^ from 
microwave conductivity measurements. The reduction of 
the real part of the conductivity ui from the Drude con- 
ductivity (Jdc and the steep upturn in the imaginary part 
of the conductivity (T2 may be interpreted in terms of 
preformed pairs leading to a pseudogap in this material. 
A pseudogap is predicted by the RVB theory of organic 
superconductivityi^i^^ 

The experimental evidence from measurements of 
1/TiT, i^s, and heat capacity all seem to point 
to the existence of a pseudogap below Tnmr in 



(ET)2Cu[N(CN)2]Br and k-(ET)2Cu(NCS)2. Thus a 
phenomenological description which takes into account 
both the spin fluctuations which are important above 
Tnmr and a pseudogap which dominates the physics be- 
low Tnmr would seem to be a reasonable starting point 
to explain the NMR data for the entire temperature 
range (clearly superconductivity must also be included 
for T < Tc). We will pursue this approach in our fu- 
ture work. In particular, if there is a pseudogap then 
important questions to answer include: (i) How big is 
the pseudogap and what symmetry does it have? (ii) 
How similar is the pseudogap in k-(ET)2X to the pseu- 
dogaps in the cuprates and in other strongly correlated 
materials such as manganites and heavy fermions? (iii) 
Is there any relationship between the pseudogap and the 
superconducting gap in k-(ET)2A? The answer to these 
questions may help put constraints on the microscopic 
theories. 

Future experiments. There are a number of key ex- 
periments required to resolve the issue whether or not 
a pseudogap is present in the low temperature metallic 
phase of k-(ET)2A". The pressure and magnetic field de- 
pendence of the nuclear spin relaxation rate and Knight 
shift will be valuable in determining the pseudogap phase 
boundary, estimating the order of magnitude of the pseu- 
dogap, and addressing the issue how the pseudogap is 
related to superconductivity. In the cuprates, there 
have been several investigations of the magnetic field de- 
pendence of the pseudogap seen in NMR experiments. 
For Bi2Sri.6Lao.4Cu06 the nuclear spin relaxation rate 
does not change with field up to 43 T.^"* However, since 
T* ~ 200 K, one may require a larger field to reduce the 
pseudogap. Similar results were found in YBa2Cu408-— 
However, in YBa2Cu307_5 [see especially Fig. 6 of Ref. 
,66j| a field of order 10 T is enough to start to close the 
pseudogap. 

The interlayer magnetoresistance of the cuprates has 
proven to be a sensitive probe of the pseudogap j^^i^^'^^i^° 
Moreover, it is found that for the field parallel to the 
layers (which means that Zeeman effects will dominate 
orbital magnetoresistance effects) the pseudogap is closed 
at a field given by 

Hpo - ^ (17) 

where T* is the pseudogap temperature. For the hole 
doped cuprates this field is of the order 100 T. In con- 
trast, for the electron-doped cuprates this field is of the 
order 30 T (and T* - 30 - 40 K), and so this is much 
more experimentally accessible. The field and temper- 
ature dependence of the interlayer resistance for several 
superconducting organic charge transfer salts^^ is qual- 
itatively similar to that for the cuprates. In particular, 
for temperatures less than the zero-field transition tem- 
perature and fields larger than the upper critical field, 
negative magnetoresistance is observed for fields perpen- 
dicular to the layers. A possible explanation is that, as 
in the cuprates, there is a suppression of the density of 
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states near the Fermi energy, and the associated pseudo- 
gap decreases with increasing magnetic field. 

A Nernst experiment can be used to probe whether 
there are superconducting fluctuations in the pseudogap 
phase, as has been done in the cuprates,-2i This experi- 
ment is particularly important in understanding the re- 
lation between the pseudogap and superconductivity. 

One could also study the pressure dependence of the 
linear coefficient of heat capacity 7. Since 7 is propor- 
tional to the density of states at the Fermi energy, a 
detailed mapping of 7(P) would be an important probe 
for studying the pseudogap. Finally, measurements of 
the Hall effect have also led to important insights into 
the pseudogap of the cuprates^^ and so perhaps the time 
is ripe to revisit these experiments in the organic charge 
transfer salts. 



V. CONCLUSIONS 

We have applied a spin fluctuation model to study the 
temperature dependences of the nuclear spin relaxation 
rate, Knight shift, and Korringa ratio in the metallic 
phase of several quasi two-dimensional organic charge 
transfer salts. This model was based on Moriya's self 
consistent renormalization theory^^ which was then ap- 
plied by Millis, Monien, and Pines^^ to cuprates. The 
large enhancement of 1/TiT between Tnmr {~ 50 K in 
K-(ET)2Cu[N(CN)2]Br} and room temperature has been 
shown to be the result of strong antiferromagnetic spin 
fluctuations. 

The antiferromagnetic correlation length is esti- 
mated to be 2.8 ± 1.8 lattice spacings in k- 
(ET)2Cu[N(CN)2]Br at T = 50 K. This value falls be- 
tween those for the Heisenberg model on the isotropic 
triangular lattice and the square lattice. 



The spin fluctuations in k-(ET)2Cu[N(CN)2]C1, k- 
(ET)2Cu[N(CN)2]Br, K(d8)-(ET)2Cu[N(CN)2]Br, and k- 
(ET)2Cu(NCS)2 are found to be similar both qualita- 
tively and quantitatively. Strong spin fluctuations seem 
to be manifested in materials close to Mott transition. 
Recent NMR experiments!?, on K-(ET)2Ag(CN)2-H20, 
which is situated further away from the Mott transition, 
suggests that the spin fluctuations in this materials are 
not as strong as those in the other k salts studied here. 

The temperature dependence of 1/TiT for T>Tnmr 
from the spin fluctuation model is qualitatively simi- 
lar with the predictions of dynamical mean fleld the- 
ory (DMFT). Below Tnmr, the spin fluctuation model 
predicts a monotically increasing 1/TiT with decreasing 
temperature while DMFT produces a constant l/TiT. 
Neither of these models can account for the large suppres- 
sion of l/Ti, Kg, and /C below Tnmr '--^ 50 K observed 
in all the K-salts studied here. This suggests two things. 
First, the low temperature regime is more complicated 
than the rcnormalized Fermi liquid previously thought to 
be the correct description of the low temperature metal- 
lic phase in these materials. Second, a pseudogap exists 
at low temperatures near the Mott insulating phase of 
the organic charge transfer salts. 
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